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Background: The Visayan Tarictic Hornbill {Penelopides panini) and the Walden's Hornbill {Aceros waldeni) are two 
threatened hornbill species endemic to the western islands of the Visayas that constitute - between Luzon and 
Mindanao - the central island group of the Philippine archipelago. In order to evaluate their genetic diversity and to 
support efforts towards their conservation, we analyzed genetic variation in ~ 600 base pairs (bp) of the 
mitochondrial control region I and at 12-19 nuclear microsatellite loci. The sampling covered extant populations, 
still occurring only on two islands {P. panini: Panay and Negros, A waldeni: only Panay), and it was augmented with 
museum specimens of extinct populations from neighboring islands. For comparison, their less endangered 
(= more abundant) sister taxa, the Luzon Tarictic Hornbill (P. manillae) from the Luzon and Polillo Islands and the 
Writhed Hornbill (A leucocephalus) from Mindanao Island, were also included in the study. We reconstructed the 
population history of the two Penelopides species and assessed the genetic population structure of the remaining 
wild populations in all four species. 

Results: Mitochondrial and nuclear data concordantly show a clear genetic separation according to the island of 
origin in both Penelopides species, but also unravel sporadic over-water movements between islands. We found 
evidence that deforestation in the last century influenced these migratory events. Both classes of markers and the 
comparison to museum specimens reveal a genetic diversity loss in both Visayan hornbill species, P. panini and 
A waldeni, as compared to their more abundant relatives. This might have been caused by local extinction of 
genetically differentiated populations together with the dramatic decline in the abundance of the extant 
populations. 

Conclusions: We demonstrated a loss in genetic diversity of P. panini and A waldeni as compared to their sister 
taxa P. manillae and A. leucocephalus. Because of the low potential for gene flow and population exchange across 
islands, saving of the remaining birds of almost extinct local populations - be it in the wild or in captivity - is 
particularly important to preserve the species' genetic potential. 
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Background 

The Philippine archipelago consists of more than 7,000 
islands, mostly of volcanic origin, with a well known 
geological and plate tectonic history [1]. According to 
various authors e.g. [2-6], during the Pleistocene glacia- 
tions, the Philippine islands were less fragmented than 
today; groups of islands formed composite Pleistocene 
Aggregate Island Complexes (PAICs) [5], that were iso- 
lated from one another by deep-water channels. Apart 
from Greater Palawan, which, biogeographically, does 
not belong to the oceanic region of the Philippines [4], 
there were three major PAICs: Greater Luzon, Greater 
Negros-Panay, and Greater Mindanao [3]. The geological 
history has had a great impact on the evolution of the 
Philippines' fauna and flora [7] resulting in an extremely 
high rate of over 57% endemism in the major faunal 
groups [8], with many lineages limited to single islands 
or PAICs [4]. Among Philippine vertebrates, this is espe- 
cially true for reptiles, amphibians, and mammals. In 
contrast, birds show, due to their volant nature, rela- 
tively low levels of endemism. Here, dispersal may be 
thus more an important driver of species diversity than 
vicariance [9]. 

Because of their species richness, the Philippines be- 
long to the 25 biodiversity hotspots on earth, but are 
simultaneously the one with the lowest percentage (3%) 
of remaining primary vegetation [10]. According to sat- 
ellite data from 1987, forest cover has declined to 23.8%, 
but varies considerably across the archipelago [11]. For 
example, Luzon retained 24% forest cover and Mindanao 
29%; in the Western Visayas, Panay retained 8% and 
Negros 4%. Even these low numbers are probably 



overestimates, as only a proportion comprises closed- 
canopy forest [11,12]. In 2002, estimates for the whole 
archipelago of 21.7% and 24.4% (depending on source) 
forest cover reflect some slowing down in the rate of de- 
forestation [13,14]. The remaining forest is of critical im- 
portance for the survival of many threatened species, 
including the hole-nesting hornbills (Bucerotidae), 
which, due to their size, depend on primary or old sec- 
ondary growth for their reproduction. 

The distribution of hornbills in the Philippines largely 
corresponds to the composition of the PAICs. These 
birds are a prominent faunal element of the Archipelago 
such that, based on the occurrence of the different horn- 
bill species, Kemp [15] divided the Philippines into six 
zoogeographical regions, the Luzonian, the Mindoroan, 
the Panay-Negrosian, the Samarian the Mindanaoean, 
and the Suluan, corresponding to the three major and 
two smaller PAICs (according to Heaney [2], the Samar- 
ian and the Mindanaoean formed one PAIC). 

The focal species of this study, the Visayan Tarictic 
Hornbill (Penelopides panini) and the Walden's Hornbill 
(Aceros waldeni), occur in the Panay-Negrosian region 
(Greater Negros-Panay, Figure 1), the most threatened 
faunal region of the Philippines [16]. For comparison, 
their less endangered (= more abundant) respective sis- 
ter taxa, the Luzon Tarictic Hornbill (P. manillae) from 
the Luzonian (Greater Luzon) and the Writhed Hornbill 
(A. leucocephalus) from the Mindanaoean (Greater Min- 
danao) are included in this study [15]. 

A. waldeni and P. panini are endemic to the western 
islands of the Visayas, which constitute - between Luzon 
and Mindanao - the central island group of the 
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Figure 1 Former and recent distribution of (A) Penelopides manillae and P. panini and of (B) Aceros waldeni and A. leucocephalus. The 

dashed lines approximately delineate the distribution ranges of the respective sister species. The illustration of the Philippines is restricted to 
larger islands and islands with former or recent occurrence of the studied hornbill species. Labeled islands are mentioned in the text. Sampled 
populations are indicated by island names in bold face. 
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Philippine archipelago. Due to deforestation and hunt- 
ing, their populations decreased rapidly [12]. P. panini 
is known from Panay (including offshore islands), 
Guimaras, Negros, Masbate, and Ticao Island (here with 
a second recognized subspecies P. p. ticaensis [17,18]). 
Today, it is still fairly common on Panay, less so on 
Negros, and extinct, or nearly so, on all other islands 
(Figure 1A) [19-21], EC pers. obs. Recently, a remnant 
population of three birds was found in Milagros, Mas- 
bate, but this very small native forest patch is slowly 
being converted into agricultural land, rendering the sur- 
vival probability of this remnant population very low 
[22]. P. panini is listed as endangered by BirdLife Inter- 
national [12]. 

A. waldeni used to occur on three islands, Panay, Gui- 
maras, and Negros (Figure IB) [18], but has survived 
with a breeding population of substantial size only on 
Panay, in the forest of the Central Panay Mountain 
Range [23]. It is extinct on Guimaras and may be extinct 
on Negros as well, as no recent sightings have been 
reported [20], EC pers. obs. It is listed as critically 
endangered by BirdLife International [12]. 

Captive breeding of both threatened Visayan hornbill 
species is very rare and incipiently. The world's only cap- 
tive population of A. waldeni (22 individuals) is held in a 
breeding facility on Panay. In 2005, for the first time, a 
total of three offspring were successfully bred by two 
pairs [24]. P. panini have bred there successfully since 
2003 [24]. Former captive reproduction success within 
the genus Penelopides could have also include P. panini, 
but due to uncertainty about species level in this genus 
[15,18], the different species most probably had been 
hybridized in captivity [25]. 

The less endangered P. manillae is classified into two 
subspecies, separated morphologically by size, the width of 
the tail band and the glossiness of the upperparts 
[15,18,26]; P. m. manillae occurs on Luzon and smaller ad- 
jacent islands, such as Marinduque and Catanduanes, P. m. 
subniger on Polillo and Patnanungan Islands (Figure 1A) 
[18,26-28]. BirdLife International [12] lists the species in 
the category 'least concern'. 



The less endangered A. leucocephalus occurs on Min- 
danao and its satellite islands (Figure IB) [18]. Its con- 
servation status is poorly known, but it is considered 
locally fairly common in suitable habitat and is listed as 
near threatened by BirdLife International [12]. 

In order to support conservation efforts for Philippine 
hornbill species, an assessment of the genetic popula- 
tion structure of the remaining wild populations is a 
basic requirement. Not many population genetic studies 
have so far been carried out in the Philippines; this is 
especially true for birds. Patterns of colonization and 
diversification in Passeriformes have been studied in a 
phylogeographic context by Jones and Kennedy [7,29] 
and Silva-Iturriza et al. [30], and at higher taxonomic 
levels by Oliveros and Moyle [9] and Sanchez-Gonzalez 
and Moyle [31]. Here, we reconstruct the phylogeo- 
graphic history of the two Penelopides species in order 
to relate it to the geological history of the Philippines. 
We also give estimates on the genetic diversity both in 
the endangered species (A. waldeni and P. panini) and 
their less threatened sister taxa (A. leucocephalus and 
P. manillae) in order to test, whether the decline in 
population sizes in the endangered species is also 
reflected in lower genetic diversity. 

Methods 

Sampling and DNA extraction 

One drop of blood was taken from captive hornbills and 
stored in Queen's Lysis Buffer [32]. The number of 
sampled individuals per population, their origin and ab- 
breviation used in the following text and figures are 
given in Table 1. In addition, 15 feather samples of P. 
panini and 7 feather samples of A. waldeni were col- 
lected from molted feathers in nest holes of wild birds 
on Panay. Furthermore, 24 museum samples from The 
Natural History Museum, Tring (U.K.; BMNH) and the 
Museum fur Naturkunde, Berlin (Germany; ZMB) - all 
dating back to the 19 th century - were included in the 
study (Table 2). 

DNA extraction of blood samples was performed 
using the DNeasy Tissue Kit (Qiagen) according to the 



Table 1 Blood samples of captive hornbills included in the study 



Number of 
samples 


Species 


Origin 


Abbreviation 
in text 


62 


P. panini 


Panay (confiscated from poaching — > of wild origin) 


PpPa 


14 


P. panini 


Negros (Bacolod City Biodiversity Conservation Center of the Negros Forest and 
Ecological Foundation, Inc/Ecological Park, Dumaguete) 


PpNe 


17 


P. manillae manillae 


Luzon* (Birds International, Inc., Quezon City/Malagos Garden Resort, 
Davao City/Avilon Montalban Zoological Park, Rizal) 


Pmm 


9 


P. manillae subniger 


Polillo* (Birds International, Inc., Quezon City) 


Pms 


22 


A. waldeni 


Panay (confiscated from poaching — » of wild origin) 


Aw 


4 


A. leucocephalus 


Mindanao (Malagos Garden Resort, Davao City/Avilon Montalban Zoological Park, Rizal) 


Al 



* origin concluded from morphological determination of captive individuals. 
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Table 2 Museum specimens included in the study. BMNH, The Natural History Museum, Tring (U.K.); ZMB, Museum fur 
Naturkunde, Berlin (Germany) 



\/r*i irh^r 

VUULI ICI 




Species 




IdanH Inralitv 
ijiaiiu, imaiiiy 


Collecting date 


vJtrl 1 Da 1 1 r\ dLL, 1 ID. 


RMNH 18Qn4 1 S 67 


Ppnp r\ninpa 
rci iciuuiucj 


nnnini 

UUI III II 


nnn in i 

UUI III II 


(~-,\ i i ma rac 
vjUI 1 1 Id 1 Qj 


f. Ian 1 999 
u Jell I I ooo 




JX273923 


RMNH 1 9Qf, A 1 ^ fSK 
Diviivn i oyo.^t. i j.oo 


Ppnpir\ninp c 
rCt iCIUpiUCi 


nn n i rt t 
jJU! Ill II 


nn mm 
put III II 


vjU 1 1 1 Id 1 di 


9R Rpr 1 9^ 

ZO UcL I Oi 




IY973Q94 
J AZ / DyZH 


RMNH 19,Qfi> 4. 1 ^ fiQ 

D I VI I N I I I OyO.H. I J. Or? 


Ppnplr\nir\pz 
rCt ICIUpiUCi 


nn n in i 
fJUl III II 


nn ntm 
put III II 


vjU 1 1 1 Id 1 di 


93 Ian 1 99 
ZD Jd I l I oo 




IY973Q9 c ; 
J AZ / DyZD 


RMNH IRRSm^DQQ 
□iv II vi i i ooo. i u. j\j.yy 


ppnp r\ninpa 
rci iciuuiucj 


nnnini 

UUI III II 


nn n in i 

UUI III II 


(~-,\ i i ma rac 
sjUI 1 1 io 1 Oj 


Mar 1879 

i vid i i o / z 




IY973Q96 

J AZ / JjZU 


RMNH 19,9,9, 10 30100 

□IV II VI 1 1 OOO. 1 \J.J)\J. 1 WW 


Ppnp r\ninpa 
rci iciuuiucj 


nnnini 

UUI III II 


nn n in i 

UUI III II 


(~-,\ i i ma rac 
UUI 1 I lal Qj 


Mar 1R79 
i vid i i o / z 




]y973Q97-iy973Q9R 
jaz/ J>y z. / jaz; J7ZO 


RMNH 1 999 1 0 3D 101 

□I VII VI 1 1 OOO. 1 U.J W. 1 W 1 


Ppnp r\r>irlpa 
rci iciuuiucj 


nnnini 

UUI III II 


nn n in i 

UUI III II 


(~-,\ i i ma rac 
sjUI 1 1 lal Go 


Mar 1879 

I Via 1 1 O / Z 




JX273929 


RMNH 1 999 1 D 30 1 09 
Divinn i ooo. i u.ju. i wz 


Ppnpir\ninp c 
rci ICIUpiUCj 


nnnini 
fJUl III II 


nn ntm 
put III II 


\-i i i i ry^ a ra c 
VJU 1 1 1 Id 1 di 


Mar 1879 
I Vld 1 1 O / Z 




IY973Q30- IY973Q31 
JAZ/ j?jU JAZ/ I 


7MR 91 77^ 

Z1IVID Z 1 / / J 


Ppnpir\nir\pz 
rci ICIUpiUCj 


nnnint 
fJUl III II 


nn ntm 
put III II 


\-i i i i a ra c 
VJU 1 1 1 Id 1 di 


Mar 1879 
I Via 1 1 O / Z 




IY973Q39 
JAZ / DyDZ 


RMNH 18Q6 4 1 S 70 

UIVllNl I IO"U,4. IJ,/u 


Ppnp r\ninpa 
rci iciuuiucj 


nnnini 

UUI III II 


nn n in i 

UUI III II 


Ma^hatp 

1 Vld jUo lc 


1 Mav 1 99 
I ividy I oo 


3 


JX273933 


RMNH 18Qn4 1 S 71 

□IVIIVI I IOjU.t. I J./ I 


Ppnplnn inp c 
rci iciuuiucj 


nnnini 

UUI III II 


nn n in i 

UUI III II 


Ma^hatp 
i vid jUd ic 


1 Ma\/ 1 99 
i i via y i oo 


3 


JX273934 


RMNH 1 9Qf, Al^ll 


Ppnpir\nir\pz 
rci ICIUpiUCj 


nnnint 
fJUl III II 


nn ntm 
put III II 


K/lachatP 
IVldbUd Lfcr 


1 Mav 1 99 
I IVldy I oo 


3 


IY973Q3^ 
JAZ / J7JJ 


7MR Qn41 

ZJV1D y\DH \ 


Ppnp\r\ninp c 
rci ICIUpiUCj 


mnn il inp 

i nut imuc 


mnn iitnp 
1 1 IUI IlllUc 


i \~7c\n fVyl^rtil^ 
LUZUI 1, IVIdl Hid 


1 R30-1 S39 

I OjU I OjZ 




IY973Q^4 

JAZ / j7Jt 


7MR 71 77? 

Z-IV1U z i / / z 


Ppnplnn inp c 
re/ iciuuiucj 


mnnil np 

1 1 IUI IIIIUC 


mnn i np 

1 1 IUI IIIIUC 


1 i i7r\n ami inn Mani a 
LUZUI 1, alUUI IU Ivlal Hid 


Ian 1R79 
Ja il io/z 




JX273955 


7MR 91 773 

Z-IV1U Z 1 / / J 


Ppnp r\ninpa 
rci iciuuiucj 


a ffin is 


a ffin is 


1 i i7r~in ami inn Mani a"* 
LUZUI 1, a 1 UUI IU Ivlal Hid 


Feb 1872 




JX273975 


7MR 91 774 

Z_l V1U Zl / / T 


ppnp Inninpz 
rci iciuuiucj 


mnnil np 

1 1 IUI IIIIUC 


mnn i np 

1 1 IUI IIIIUC 


1 i i7r~\n ami inrl Mani a 
LUZUI 1, alUUI IU Ivlal Hid 


Jan 1872 




]Y973QSfS- IY973QS7 

JAZ/ J7JU JAZ/ J>z?J/ 


ZMB 32616 


Penelopides 


man Mae 


maniliae 


Culion, Palawan? (see discussion) 


? 




JX273958-JX273959 


ZMB 2000.22207 


Penebpides 


moniilae 


man Mae 


Luzon, Manila 


4-9 Jun 18 


\6 


JX273960 


ZMB 2000.22208 


Penelopides 


man Mae 


man Mae 


Culion, Palawan? (see discussion) 


? 




JX273961-JX273962 


BMNH 1896.4.15.98 


Aceros 


waideni 




Guimaras 


28 Jan 188 


8 


JX273808 


BMNH 1896.4.15.99 


Aceros 


waideni 




Guimaras 


28 Dec 18£ 


M 


JX273809 


BMNH 1896.4.15.100 


Aceros 


waideni 




Guimaras 


27 Dec 1 8f 


M 


JX273810 


ZMB 2000.22212 


Aceros 


waideni 




Guimaras 


6 Jan 1888 




JX27381 1 


ZMB 2000.22215 


Aceros 


waideni 




Guimaras 


Jan 1888 




JX273812 


BMNH 1897.5.13.471 


Aceros 


waideni 




Negros 


27 Mar 1& 


6 


JX273813 



* for a note on taxonomic status and origin stated on the museum label of ZMB 21 773 see Additional file 1 . 



manufacturer's instructions for blood samples. Feather 
samples were digested with 15 ul proteinase K (20mg/ 
ml) in the presence of 40 ul DTT and 5 ul Carrier 
RNA (1 ug/ul) in 450 ul buffer solution (100 mM Tris, 
lOmM EDTA, lOOmM NaCl, 0.1% SDS). After incuba- 
tion at 55°C over night, we added again 15 ul protein- 
ase K (20mg/ml) and incubated the mix for one more 
hour. Afterwards, we purified the samples according to 
the manufacturer's instructions of the DNeasy Tissue 
Kit (Qiagen), but increased EtOH to 400 ul and the 
washing solutions to 600 ul. DNA was eluted twice in 
50 ul of the buffer solution provided with the kit. 

DNA extraction from toe-pad samples of museum 
specimens were performed using the DNeasy Tissue 
Kit (Qiagen), according to the manufacturer's instruc- 
tions, adding 5 ul Carrier-RNA (1 ug/ul). In all steps, 
solution volumes were doubled. DNA was eluted twice, 
initially in 50 ul and then in 30 ul of the kit buffer solu- 
tion. DNA extraction from museum samples was car- 
ried out in a dedicated ancient DNA laboratory, 
physically separated from the laboratory where contem- 
porary DNA was handled. 



mtDNA amplifying and sequencing 

The mitochondrial genome of Philippine hornbills con- 
tains two control regions of which we amplified the 5'end 
of the first control region (CRI) as described for fragment 
3 in Sammler et al. [33]. Products were sequenced with 
primer AcePen_Glu-for to obtain the targeted first ~ 650 
base pairs. 

Since the amount of DNA isolated from molted feathers 
and museum specimens is low and the DNA often 
degraded, the corresponding sequences of these samples 
were amplified in two short, partially overlapping frag- 
ments. The first fragment was amplified with the forward 
primer AcePen_Glu-for and a CRI specific reverse primer 
(Penpan_CRI-spec-180rev or PenpanGui_CRI-spec-182rev 
for P. panini, Penmanman_CRI-spec-185rev for P. m. 
maniliae, Penaff_CRI-spec-186rev for museum specimen 
ZMB 21773, and Acewal_CRI-spec-180rev for A. wai- 
deni). The second fragment was amplified with a CRI 
specific forward primer (Penpan_CRTspec-128for for P. 
panini, Penmanman_CRTspec-115for for P. m. maniliae, 
Penaff_CRI-spec-124for for museum specimen ZMB 
21773, and Acewal_CRI-spec-130for for A, waideni) 
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and the reverse primer AcePen_644rev or, when 
amplificates did not yield sufficient results, also with 
AcePen_426rev (Table 3). 

The 37.5 ul PCR reaction volumes (buffer solution: 
lOmM Tris-HCl, pH 9.0, 50mM KC1, 1.5mM MgCl 2 , 
0.1% Triton X100, 0.2 mg/ml Bovine Serum Albumin 
(BSA)) were set up as follows: 6 ul dNTP Mixture (2.0 
mM each), 1.2 ul of each primer (10 uM), 3 ul DNA tem- 
plate, 0.25 ul Taq polymerase (5U/ul, MP Biomedicals). 
The reaction was performed under the following condi- 
tions: initial denaturation at 94°C for 5 min; 50 cycles: 
94°C for 90 s, 50°C for 1 min, 72"C for 1 min; and a final 
extension at 72°C for 10 min. Products were sequenced 
as described in Sammler et al. [33] with AcePen_Glu-for, 
AcePen_644rev or AcePen_426rev, respectively. 

To verify the authenticity of amplificates from histor- 
ical samples, all the mtDNA PCRs of museum samples 
were repeated by another person in a separate laboratory 
in a plant genetics department, where animal DNA had 
never been handled. Specifically, all repeated PCR/ 
sequencing analyses of the first fragment yielded identical 
results in 1 st and 2 nd PCR. For the second fragment, 1 st 
and 2 nd PCR also yielded identical results, except for 
two occasions with a single bp mismatch. Scarcity of 
DNA precluded a third PCR, but these two mismatches 
translate into an per base pair error rate of about 
0.014%, a measure in the range of known Taq error 
rates. There was no indication of an elevated error rate 
in the historical samples. 

Sequences were aligned in BioEdit version 7.0.5.3 [34]. 
Mitochondrial haplo types were defined based on a 591- 
593 bp data partition. Distance matrices and statistical 
parsimony haplotype networks were constructed using 
TCS 1.21 with default parameter settings, but a fixed 
connection limit at 100 steps [35]. 



Microsatellite genotyping 

All blood samples were genotyped at 12-19 previously 
published polymorphic microsatellite loci (Table 4). PCR 
and fragment size analysis of listed loci were carried out 
as described in Sammler et al. [36]. Annealing tempera- 
tures of Bbi and Bubi primers were set to 55°C. 

mtDNA data analysis 

Standard measures of genetic diversity (haplotype diver- 
sity, HD [38]; nucleotide diversity, ND [39]) as well as of 
genetic divergence among island populations (fixation 
index, F sx [40]), were calculated as implemented in the 
software package ARLEQUIN version 3.5.1.2 [41]. When 
more than two populations were compared (i.e., P. pa- 
nini), we performed an analysis of molecular variance 
(AMOVA) with the same software. Significance of F sx / 
AMOVA results was evaluated by permutation analysis 
(1,000 permutations). Diagnostic mutations of each island 
population were identified visually in the alignment. 
ARLEQUIN was also used to assess the fit of the demo- 
graphic expansion model [42] to hornbill genetic varia- 
tions (each population singularly). Mismatch distributions 
were used to fit the model of sudden demographic expan- 
sion. Goodness of fit was assessed by the sum of square 
deviations (SSD) between the observed and the expected 
mismatch and its significance (P-SSD) was determined by 
a parametric bootstrap with 10,000 replicates. Mismatch 
distributions were not calculated for A. leucocephalus be- 
cause of the small sample size of the single population of 
this species included in the study. 

Microsatellite data analysis 

The program FSTAT 2.9.3.2 [43] was used to calculate 
allelic richness (AR). F sx values were calculated with 
ARLEQUIN. Significance of these values was evaluated 



Table 3 PCR primers used to amplify CRI gene fragments. Y = C/T; R= A/G 



PCR primers 


Primer sequence (5'-3') 


Reference 


AcePen_Glu-for a 


GCT TIT CTC CAA GGT CTA CAG CTC 


[33] 


AcePen_Cyt1018-rev a 


GGG TGT TCT ACT GGT TGG CTG CC 


[33] 


Penpan_CRI-spec-1 80rev 


CRY YRT YYA CAT TAA GTG A 


This study 


PenpanGui_CRI-spec-182rev 


GTC ATT ATC TAC ATT AAG TAA G 


This study 


Penmanman_CRI-spec-185rev 


TTA ATA TGT CGT TGT TTA CAT G 


This study 


Penaff_CRI-spec-186rev 


GTA TGT CAT TGC TTG CAT TGA GTA G 


This study 


Acewal_CRI-spec-1 80rev 


CAT TRT CTG CAT TTA AGC GT 


This study 


Penpan_CRI-spec-1 28for 


ACG ACT AGT TAT TAA TGC T 


This study 


Penmanman_CRI-spec-1 1 5for 


CAT AAG GTA ATG CTC TAT ACA ATT 


This study 


Penaff_CRI-spec-1 24for 


CTA TAT GAT TAA CTA TTA ATG CTC 


This study 


Acewal_CRI-spec-1 30for 


GAT TGA CTG TCA ATG TTT GT 


This study 


AcePen_644rev a 


AAG GGA ACC AAC AGT GCC AAA C 


This study 


AcePen_426rev a 


GTT GCT GAT TTC TCG TGA GG 


This study 



a primers with two annealing sites due to the duplication of a large fragment of the hornbills' mt genome [33]. 
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Table 4 Microsatellite loci used in the study 



Loci 


Genotyped species 


Reference 


Bbi2 


P. panini, P. manillae, A. waldeni, 
A. leucocephalus 


[37] 


Bbi7 


P. panini, P. manillae, A. waldeni, 
A. leucocephalus 


[37] 


Bbi 1 3 


P. panini, P. manillae 


[37] 


Bbi16 


P. panini, P. manillae, A. waldeni, 
A. leucocephalus 


[37] 


Bubi294 


P. panini, P. manillae, A. waldeni, 
A. leucocephalus 


[37] 


Bubi326 


P. panini, P. manillae, A. waldeni 


[37] 


Pp_GA_3 


P. manillae 


[36] 


Pp_GA_4 


P. panini, P. manillae 


[36] 


Pp_GT_1 


P. panini, P. manillae, A. waldeni, 
A. leucocephalus 


[36] 


Pp_GT_3 


A. leucocephalus 


[36] 


Pp_GT_4 


P. panini, P. manillae, A. waldeni, 
A. leucocephalus 


[36] 


Pp_GT_5 


P. panini, P. manillae 


[36] 


Pp_GT_6 


P. panini, P. manillae, A. waldeni, 
A. leucocephalus 


[36] 


Pp_GT_7 


P. panini, P. manillae, A. waldeni, 
A. leucocephalus 


[36] 


Pp GT 10 


r. uuillill, r. nlUfllllUc, At. wuiucl 11 




Pp_GT_1 7 


P. panini, P. manillae 


[36] 


Pp_GT_1 8 


P. panini, P. manillae, A. waldeni, 
A. leucocephalus 


[36] 


Pp_GT_1 9 


P. panini, P. manillae, A. waldeni, 
A. leucocephalus 


[36] 


Pp_GT_21 


P. panini, P. manillae, A. waldeni, 
A. leucocephalus 


[36] 


Pp_GT_22 


P. panini, P. manillae 


[36] 



by permutation analysis (1,000 permutations). Private 
alleles (PA) were assessed with the Excel Add-In soft- 
ware GENALEX 6 [44]. ARLEQUIN was used to calcu- 
late observed and expected heterozygosity (H Q and H E , 
respectively) and to test for deviation from Hardy- 
Weinberg equilibrium (HWE) for each locus in each 
population using Fisher's exact test and the Markov 
chain method (1,000 demorization steps, 100 batches, 
with 10,000 iterations per batch set). Levels of signifi- 
cance of the HWE and linkage disequilibrium tests were 
Bonferroni corrected for multiple comparisons. A Bayes- 
ian clustering was performed using the software 
STRUCTURE version 2.3.2 [45]. Genetic subdivision was 
evaluated estimating the likelihood and sample compos- 
ition of independent runs of subgroups (K = 1-5), as- 
suming an admixture model, with a burn-in length of 
100,000 and a data collection period of 900,000 itera- 
tions. To check for convergence of the Markov Chain 
parameters, ten replicate runs for each value of K were 
performed. The value of K best suited for the data was 



determined according to Evanno et al. [46] and Earl 
and von Holdt [47]. 

Demographic signatures of bottlenecks were investi- 
gated using the BOTTLENECK 1.2.02 software [48] under 
three different possible evolutionary models, i.e., Infinite 
Allele Model (IAM), Stepwise Mutation Model (SMM), 
and Two-Phase Model (TPM). Estimations of significant 
deviations from the null hypothesis of mutation-drift 
equilibrium were based on 10,000 replications. 

Results 

mtDNA diversity and divergence 
Penelopides panini 

Based on the first 593 bp of the mitochondrial control re- 
gion I sequences obtained for 102 specimens, we found 
76 polymorphic sites (one indel, 75 transitions), defining 
62 unique haplotypes (Figure 2A) [GenBank: HQ834451, 
HQ834457-HQ834461, JX273818-JX273935]. Due to sin- 
gle heteroplasmic sites cf. [33,49], 22 specimens exhibited 
two haplotypes. The frequency of such haplotypes was 
weighted by 0.5 in the frequency representation of Fig- 
ure 2. The two most abundant haplotypes were found in 
8 individuals each; all others are only represented by four 
or less individuals (Figure 2A). The AMOVA revealed a 
highly significant overall divergence among populations 
(F ST =0.572, p<0.001), such that 57% of the mtDNA vari- 
ation was apportioned to divergence among populations. 
All single populations had significantly diverged from 
each others (pairwise F ST between 0.284 and 0.664; 
p<0.01; see below for details). Three individuals originat- 
ing from Negros (PpNell, 40, 55) carry two unique hap- 
lotypes that differ by 18-24 substitutions from all other 
individuals coming from the same area. Both haplotypes 
are genetically more similar to haplotypes originating 
from Panay/Guimaras (4 and 5 substitutions away from 
the respective closest haplotype, Figure 2A). Apart from 
these three individuals, the two recent populations of 
Panay and Negros are clearly separated by at least 17 mu- 
tational steps in the network and seven diagnostic muta- 
tions in the alignment. The fixation coefficient (F S t) 
between these populations is 0.665 (p<0.001), but rises to 
0.765 (p<0.001) when calculated excluding the three indi- 
viduals exhibiting deviant haplotypes (Ppll/Pp40 and 
Pp55). Haplotype diversities range from 0.900 (Negros) to 
1.000 (Masbate). The Negros population has a nucleotide 
diversity higher than that found on Panay (ND = 0.015 
vs. 0.010, Table 5), but when the two exceptionally diver- 
gent haplotypes (PpNell/40, PpNe55) are removed from 
computations, ND drops to 0.004. The highest number 
of base substitutions between two individuals of the 
species is 28, corresponding to a nucleotide divergence 
of 4.7%. 

The museum specimens from Guimaras and Masbate 
are all represented by unique haplotypes that cluster closer 
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A- Penelopides panini 




B- Penelopides manillae 
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Figure 2 Haplotype networks of (A) Penelopides panini, (B) P. manillae, (C) Aceros waldeni, and (D) A leucocephalus. For heteroplasmic 
individuals, each haplotype frequency was weighted by 0.5. Italic labels indicate individuals that are discussed in detail in the text. Haplotypes 
based on museum specimens are marked by an asterisk, those of P. manillae with known origin from the Manila area additionally by an M. 



to haplotypes from Panay than to those from Negros 
(Figure 2A). There are no diagnostic mutations in the align- 
ment which would clearly separate them from the Panay 
population. However, individuals of those two populations 
do not cluster randomly between the Panay haplotypes, but 
form separate subclusters (Figure 2A), 1-7 and 4-7 muta- 
tional steps away from the next Panay haplotype. Fixation 
indices to Panay are lower (Panay/Guimaras 0.295, 
p<0.001; Panay/Masbate 0.284, p<0.001) than to Negros 
(Negros/Guimaras 0.640, p<0.001; Negros/Masbate 0.576, 
p<0.001). Mismatch distributions do not deviate from the 



purely demographic expansion model (0.288 < P-SDD < 
0.773) with the only exception of the population from Neg- 
ros (P-SDD = 0.048). 

Penelopides manillae 

Based on 591 bp sequence in 32 specimens, we found 98 
polymorphic sites (1 transversion, 96 transitions, 1 site 
with both a transversion and a transition), defining 26 
unique haplotypes (Figure 2B) [GenBank: JX273936- 
JX273974]. Due to heteroplasmic sites, seven specimens 
are assigned to two haplotypes. 



Table 5 mtDNA and microsatellite diversity measures 



P. panini 
(Panay) 



P. panini 
(Negros) 



P. panini 
(Guimaras) 



P. panini 
(Masbate) 



P. manillae manillae 
(Luzon) 



P. manillae subniger 
(Polillo) 



A. waldeni 
(Panay) 



A. waldeni 
(Guimaras) 



A. waldeni 
(Negros) 



A. leucocephalus 
(Mindanao) 



mt DNA 

n/H 77/43 

HD 0.975±0.005 

ND 0.010±0.006 

DM 0 

microsat 

loci 18 

n 62 

H 0 0.586 ± 0.280 
H E 0.640 ± 0.239 
PA 2.111 
AR 5.07 



14/8 

0.900±0.046 
0.015±0.008 
T 



0.642 ± 0.236 
0.665 ±0.182 
0.611 
4.78 



8/6 

0.933±0.077 
0.008±0.005 
0 



3/3 

1 .000±0.272 
0.014+0.01 1 
0 



23/19 

0.966±0.020 
0.037±0.019 



19 
1/ 

0.547 ±0.311 
0.605 ± 0.280 
3.789 
5.20 



9/7 

0.910±0.035 
0.007±0.004 



0.576 ± 0.337 
0.578 ± 0.240 
0.737 
3.53 



29/13 

0.897±0.035 
0.011 ±0.006 
0 

13 

22 

0.496 ± 0.278 
0.543 ± 0.236 

2.88 



5/3 

0.700±0.218 
0.005±0.004 
0 



1/1 



4/4 

1.000+0.177 
0.038±0.026 



0.625 ± 0.271 
0.645 ± 0.240 

3.91 



n number of samples analyzed, H number of haplotypes, HD haplotype diversity (± standard deviation, SD); ND nucleotide diversity (± SD), DM diagnostic mutations, H 0 mean observed heterozygosity (± SD); H E mean 
expected heterozygosity (± SD), PA private alleles, AR allelic richness over all loci. 
a when PpNe55, 1 1, 40 excluded. 
b when Pmm28 excluded. 
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No haplotype is dominant in frequency (Figure 2B). 
Apart from individual Pmm28, both subspecies, P. m. 
manillae und P. m. subniger, are separated by at least 18 
base substitutions in the network and can be distin- 
guished from one another by 7 diagnostic mutations in 
the alignment. The F ST between the two subspecies is 
0.502 (0.529 when excluding Pmm28; both values are 
significant, p<0.001). P. m. subniger has a lower haplotype 
and nucleotide diversity than P. m. manillae (HD=0.910, 
ND = 0.007 vs. HD=0.966, ND=0.037, Table 5). The high- 
est number of base substitutions between two individuals 
of the species is 70 (11.8% sequence divergence). 

In the network, P. m. manillae forms three hap- 
logroups (I, II, III, Figure 2B) separated by large genetic 
gaps. The first (I) is formed by four museum specimens 
and Pmml7. Two of the museum specimens (ZMB 9641, 
ZMB 21772) come from the Manila area, the other two 
(ZMB 2000.22208, ZMB 32616) are labeled "Culion", an 
island north of Palawan. Seventeen substitutions separate 
haplogroup I from the second one (II). This is formed by 
most of the individuals sampled, including two museum 
specimens (ZMB 21774, ZMB 2000.22207), which also 
come from the Manila area. The third haplogroup (III), 
formed by Pmm24/25, is separated by at least 52 substi- 
tutions from haplogroup II. This genetic divergence is by 
far larger than that between the two morphologically 
differentiated subspecies P. m. manillae and P. m. subni- 
ger (Figure 2B). 

Similar to what was observed in P. panini, we found 
also in P. manillae one individual (Pmm28) carrying a 
haplotype that does not fit with its morphological identi- 
fication. Pmm28 is assigned morphologically to P. m. 
manillae, but its haplotype clusters closer to those of 
P. m. subniger (4-11 base substitutions), while there 
are at least 21 substitutions between Pmm28 and the 
closest P. m. manillae haplotype. Mismatch distributions 
never deviate from the purely demographic expansion 
model (0.156 < P-SDD < 0.371). 

Aceros waldeni 

Based on 593 bp sequence of 35 specimens, we found 33 
polymorphic sites (2 indels, 1 transversion, 30 transitions), 
defining 17 unique haplotypes (Figure 2C) [GenBank: 
HQ834450, HQ834452-HQ834456, JX273781-JX273813]. 
Due to heteroplasmic sites, four specimens are assigned to 
two haplotypes. 

With nine individuals, one haplotype is the most abun- 
dant; all others are only represented by four or less indivi- 
duals (Figure 2C). mtDNA diversity is lower in Guimaras 
(HD=0.700, ND=0.005) than in Panay (HD=0.897, 
ND=0.011). The highest number of base substitutions be- 
tween two individuals is 16 (2.7% sequence divergence). 
The five museum specimens originating from Guimaras 
carry three unique haplotypes that cluster among the 
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Figure 3 STRUCTURE analysis of (A) Penelopides panini over 18 
loci and of (B) P. manillae over 19 microsatellite loci. Each 
individual is represented by one vertical column divided into (A) 
K = 2 and (B) K = 3 colored clusters. Each color line is proportional 
to the individual's estimated coefficient of membership to that 
particular cluster. Labeled individuals are discussed in detail in 
the text. 



Panay haplotypes. There are no diagnostic mutations in 
the alignment which would separate the individuals 
from Panay and Guimaras; the fixation index is 0.084 
(p=0.063). The haplotype of the single museum speci- 
men originating from Negros (BMNH 1897.5.13.471) is 
separated from the next Panay haplotype by 9 muta- 
tional steps in the network (Figure 2C) and 4 diagnos- 
tic mutations in the alignment (Table 5). Mismatch 
distributions never deviate from the purely demo- 
graphic expansion model (0.393 < P-SDD < 0.373). 



Aceros leucocephalus 

Based on 592 bp sequence in 4 specimens, we found 41 
polymorphic sites (all transitions). All individuals bear 
unique haplotypes (Figure 2D) [GenBank: JX273814- 
JX273817]. HD is 1.000 and ND is 0.038, both values 
however associated with high uncertainty (as reflected in 
SD values), due to small sample size. The highest num- 
ber of base substitutions between two individuals is 30 
(5.1% sequence divergence). 

Microsatellite diversity and divergence 

The characteristics of the 14 Pp loci and their usability 
for population genetics analyses are presented by Samm- 
ler et al. [36]. Six of the eight analyzed Bbi and Bubi loci 
[37] were polymorphic, at least in the genus Penelopides. 
Their number of alleles per locus varied from 1 to 8 in A. 
waldeni, 1 to 6 in A. leucocephalus, 2 to 13 in P. panini, 
and 2 to 19 in P. manillae. We confirmed that no linkage 
disequilibrium could be found between any pair of loci. 
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Penelopides panini 

STRUCTURE results of independent runs with K ranging 
between 1 and 5 supported two as the most likely num- 
ber of clusters. The clear division into two groups 
(Figure 3A) corresponds to the two populations from 
Panay and Negros. Intriguingly, the individuals carrying 
mtDNA haplotypes deviant from their region of origin 
(PpNell/40, PpNe55; Figure 2A) are inconspicuous with 
regard to their microsatellite genotype. Interestingly, in 
two individuals from Panay (PpPal34 and PpPa78), about 
25% of their genotypes is assigned to the Negros cluster. 

The number of private alleles and the allelic richness 
over all loci is higher in the Panay than in the Negros 
population (PA = 2.11 vs. 0.61; AR = 5.07 vs. 4.78, 
Table 5) while this pattern is reversed for the observed 
and expected heterozygosities (Table 5). In the Negros 
population no loci deviate from the HWE while two loci 
are significantly out of HWE in the Panay population 
due to a lack of heterozygotes (data not shown). The fix- 
ation index between these two populations is 0.085 
(p<0.001). For the Panay population, BOTTLENECK 
found a normal L-shaped allele frequency distribution as 
expected under mutation-drift equilibrium, suggesting 
that the effective population size has remained constant 
in the past [48]. For the Negros population, however, a 
shifted mode is indicative of a recent bottleneck. 

Penelopides manillae 

STRUCTURE was run with K = 1-5 and supports three 
as the most likely number of genetic clusters in this spe- 
cies. Individuals morphologically determined as P. m. 
subniger build a first cluster, and individuals of P. m. 
manillae build a second cluster, except Pmm24 and 
Pmm25, which group on their own (Figure 3B). Pmm24 
and Pmm25 also had extremely divergent mitochondrial 
haplotypes (Figures 2B and 3). Four individuals (Pmsl, 
Pms9, Pms6, Pms3) can be assigned partially both to the 
subniger- and the manillae-cluster. The individuals 
Pmml7 (haplogroup I of P. m. manillae; Figure 2B) and 
Pmm28, clustering closer to P. m. subniger in the 
mtDNA network (Figure 2B), do not stand out at the 
microsatellite genotypes, but are genetically similar to 
the other individuals assigned to P. m. manillae on mor- 
phological grounds. 

The number of private alleles and the allelic richness 
over all loci is higher in P. m. manillae than in P. m. 
subniger (PA = 3.789 vs. 0.737; AR = 5.20 vs. 3.53, 
Table 5). H Q is nearly the same in the two subspecies 
while H E is slightly higher in P. m. manillae than in 
P. m. subniger (Table 5); one locus deviates from HWE in 
either subspecies due to heterozygote deficiency (data 
not shown). The fixation index between these two popu- 
lations is 0.050 (p<0.001). For P. m. manillae, computa- 
tions with BOTTLENECK found a normal L-shaped 



allele frequency distribution as expected under mutation- 
drift equilibrium (constancy in the past population effect- 
ive size) [48], whereas for P. m. subniger a shifted mode 
indicates a recent bottleneck. 

Aceros waldeni 

Analyses with STRUCTURE did not find any subpopula- 
tion structure in the sampling on Panay (K = 1 got best 
support; data not shown). H E is slightly higher than H Q in 
the species (Table 5) and one locus shows deviations from 
HWE due to a lack of heterozygotes (data not shown). 
Computations with BOTTLENECK found a normal L- 
shaped allele frequency distribution as expected under 
mutation-drift equilibrium [48] . 

Discussion 

Penelopides panini 

Whereas the microsatellite data for P. panini show a 
clear division between the two populations from Panay 
and Negros (the sole exceptions being PpPal34 and 
PpPa78, Figure 3A), mtDNA revealed that three indivi- 
duals from Negros (PpNell, 40, 55) carry two haplo- 
types that are genetically more similar to haplotypes 
found on Panay, Guimaras and Masbate than they are to 
the remaining haplotypes found on Negros (Figure 2A). 
If we do not consider these two deviant haplotypes and 
the microsatellite genotypes of PpPal34 and PpPa78, the 
overall clear-cut subdivision in the genetic architecture 
of the species would be largely consistent with the geo- 
logical history of the area. We hypothesize that substan- 
tial gene flow between Panay/Guimaras and Negros 
ceased at the end of the last Pleistocene ice age when 
the Panay-Negrosian PAIC broke down into separated 
islands due to eustatic changes in the sea level. While 
such a Pleistocene sea level-related model of speciation 
is rejected in favor of an earlier, Pliocene diversification 
for, e.g., rodents of the genus Apomys [50] or for the 
fruit bat Haplonycteris fischeri [51], the latter author 
suggests that sea-level changes may have played an im- 
portant role for more recent divergence at the within 
PAICs level, such as among Negros and Panay. The 
shortest distance between the two islands today is about 
15 km but, due to the presence of Guimaras and several 
other intervening islets, the longest open water distance 
is just 5 km. Tarictic hornbills are not strong flyers, and 
group-living, territorial species often make poor candi- 
dates for long-distance dispersal, even to offshore islands 
[52]. However, the three birds PpNell, 40, 55 presum- 
ably descended from at least two female ancestors, 
which had recently, i.e., some generations ago, migrated 
from Panay or Guimaras to Negros. Since their haplo- 
types cluster relatively close to haplotypes from both 
islands (Figure 2A), the exact origin cannot be ascer- 
tained. Due to the general maternal inheritance of the 
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mtDNA in vertebrates including birds (although, in rare 
occasions, paternal leakage has been found [53]), the 
haplotypes of these putative immigrants had persisted in 
the Negros population (Figure 2A), while most of their 
Panay/Guimaras ancestry has vanished at the microsatel- 
lite level due to the bi-parental inheritance of these mar- 
kers (Figure 3A). 

In the last 100 years, both Panay and Negros suffered 
from severe deforestation; Guimaras became completely 
denuded [11,14]. This complete deforestation might have 
forced birds from Guimaras to abandon the island and 
to seek still extant forested areas on Negros and Panay. 
These islands would then still host a (partial) mitochon- 
drial heritage of Guimaras, although the maternal 
lineages carried by all 19 th century Guimaras museum 
specimens could not be found in our sampling. These 
Guimaras lineages have probably gone extinct. However, 
in this regard, a Guimaras origin of PpNell/40, 55 
seems the most plausible. 

Today, the shortest distance between the remaining 
forest patches of Panay and Negros is about 100 km, a 
barrier that has become increasingly difficult for forest 
species to cross. As a consequence, we assume that gene 
flow between the two areas is currently strongly reduced 
or completely interrupted. The clear subdivision in the 
microsatellite data (Figure 3A) supports this hypothesis, 
the sole exceptions being PpPal34 and PpPa78. 

Interestingly, for another forest bird species, the Phil- 
ippine bulbul Hypsipetes philippinus, despite a small 
sample size of only three individuals from Negros, 
mtDNA suggested multiple independent colonization 
events for this island and a lack of a direct link between 
populations from Panay and Negros [30]. These authors, 
however, did not analyze any nuclear marker. A similar 
situation, again based only on mtDNA, was found for 
four out of the five passerine species studied by Jones 
and Kennedy [7,29]. However, in a fifth species, the Is- 
land Verditer-Flycatcher Eumyias panayensis, the popu- 
lations of Negros and Panay are clearly separated and 
form sister clades [7]. Generally, dispersal across islands 
appears much more infrequent in our hornbill species 
than in the more vagile passerines. 

Whereas the population of P. panini on Panay shows 
no indication of a recent bottleneck, the one on Negros 
does. This comes along with lower values of genetic di- 
versity in the Negros population (Table 5). Negros 
retained half of the forest cover of Panay (4% vs. 8%) 
[11]; in addition, the forest cover on Panay (located 
mainly in the Central Panay Mountain Range) is almost 
continuous, while it is quite fragmented on Negros. We 
are aware that our sample size for this latter population 
is rather small for a precise estimate of its genetic diver- 
sity. We do not have access to 30 individuals, which is 
recommended for bottleneck analyses by Luikart et al. 



[54]. Small sample sizes are likely to miss alleles at low 
frequency [55] and might cause allele frequency distribu- 
tions to resemble those of a bottlenecked population. 
Nonetheless, both our genetic analysis and the ecological 
factors are consistent with a declining and genetically 
depauperate population on Negros. 

Interestingly, not only the historic museum specimens 
from Guimaras, but also the three from Masbate cluster 
close to the Panay haplotypes (Figure 2A). The Negros 
population is thus the only one clearly separated from 
the other insular populations (except for three indivi- 
duals: PpNell, 40, 57). Unfortunately, we could not 
genotype the museum specimens at the microsatellite 
loci; mtDNA, however, suggests genetic connectedness 
among Guimaras, Panay and Masbate. In the case of 
Guimaras, this is not surprising; only a narrow water 
channel of about 2 km separates nowadays this island 
from Panay (Figure 1). Obviously, this channel had not 
acted as such a strong geographic barrier as the sea be- 
tween Panay/Guimaras and Negros. In the case of Mas- 
bate, it is noteworthy that the over water distance 
between Panay and Masbate is at least 35 km and thus 
longer than that between Negros and Panay. We 
hypothesize that the small islets lying in a row between 
Panay and Masbate could have acted as stepping stones 
facilitating dispersal. 

Similarly to what was found for the Guimaras speci- 
mens, the Masbate specimens also carry unique haplo- 
types. Thus, Masbate also hosts an independent (yet 
closely related to Panay) lineage, which still survives with 
a few last individuals [19,22] or might have already got 
lost, if the population went extinct since the last census. 

Penelopides manillae 

Both mtDNA and microsatellites indicate a separation 
between the two recognized subspecies P. m. manillae 
and P. m. subniger. However, divergence within one sub- 
species (manillae) is comparable to divergence between 
them (see Figure 2B and Table 5). In the haplotype net- 
work (Figure 2B), Pmm28 is an outlier. Although mor- 
phologically assigned to P. m. manillae, it clusters within 
the P. m. subniger haplogroup. The microsatellite profile 
of Pmm28 (coefficient of membership to manillae of 
about 90%; Figure 3B) indicates that this bird descends 
from a female ancestor that migrated from Polillo to 
Luzon generations ago. 

Haplogroup I in the network (Figure 2B) is formed by 
Pmml7 and four museum specimens. Two of these mu- 
seum specimens (ZMB 2000.22208, ZMB 32616) are la- 
beled "Culion", an island north of Palawan. For lack of 
any historical evidence for an occurrence of P. manillae 
on Palawan and its satellite islands, which, biogeographi- 
cally, even do not belong to the Philippines, but to 
Borneo [4], we cannot rule out the possibility of these 
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specimens having been transferred to Culion by humans. 
The correct geographic origin of these two individuals is 
thus considered uncertain. The other two (ZMB 9641, 
ZMB 21772) come from or around Manila. Two further 
museum specimens coming also from the Manila area 
(ZMB 21774, ZMB 2000.22207) cluster in haplogroup II. 
This implies that in the Manila area two divergent mito- 
chondrial lineages coexist. Such a subdivision is not evi- 
dent at the microsatellite level, though (Figure 3B). We 
cannot rule out the possibility that these results might 
be biased by our relatively small sample size that could 
have resulted in missing potentially extant intermediate 
haplotypes. On the other hand, given the strikingly 
contrasting pattern yielded by microsatellites, we tend 
to favor a scenario where two distinct mitochondrial 
lineages have coexisted for an extended period of time, 
either because of ancient larger population size or ori- 
ginating from two formerly separated populations 
within Luzon which later met and merged in the area 
around Manila. 

In contrast, haplogroup III (Pmm24/25) stands out in 
both the mitochondrial and nuclear data sets, suggesting 
the occurrence of a further distinct genetic lineage on 
Luzon. P. manillae has a very wide distribution. Although 
encompassing islands such as Marinduque and Catan- 
duanes which are as far away from the main island of 
Luzon as the Polillo Islands are, no further subspecies for 
those islands have been identified morphologically. Our 
data show that the lineage represented in our sampling by 
Pmm24/25 is genetically differentiated from both P. m. 
manillae (to which it was assigned in the first place) and 
P. m. subniger. More importantly, it is as distant from 
P. m. manillae as the latter subspecies is from P. m. 
subniger. While the instances of deviant mt haplotypes 
with inconspicuous microsatellite genotypes are indica- 
tive of historical gene flow (see above), the microsatel- 
lite genotypic compositions of Pmsl, Pms9, Pms6, 
Pms3 might suggest ongoing gene flow, as these spe- 
cimens presumably represent direct hybrids between a 
P. m. subniger female (inferred from the mtDNA type) 
and a P. m. manillae male. Although these samples 
originate from a non-breeding facility, we cannot ex- 
clude with certainty that these hybrids have been pro- 
duced in captivity [25], such that our inference about 
naturally occurring gene flow and hybridization should 
be considered preliminary. 

Apart from these hybrids, we hypothesize that P. m. 
manillae and P. m. subniger became separated when the 
sea level rose at the end of the Pleistocene and the 
Polillo Islands became disconnected from the main is- 
land of Luzon. While these islands were land-bridge 
connected to southern rather than central Luzon [2], 
they are nowadays about 20 km away from central 
Luzon. In spite of deforestation on Luzon, the coastline 



facing the Polillo Islands is still afforested, and migration 
between the islands would be still possible. Although 
tarictic hornbills are not renowned for long-distance dis- 
persal, even to nearby offshore islands [52], migration 
events become the more likely, the more the pressure on 
the Polillo and Patnanungan populations increases due 
to logging, mining, and conversion of primary forest into 
coconut plantations [28]. Thus, gene flow might have 
been directed from the small Polillo Islands to the large 
island of Luzon by individuals attempting to escape the 
increasing human pressure. Such a scenario would ex- 
plain the observed phylogeographic pattern; this, in turn, 
would be at odds with the theory of island biogeography 
that presumes large islands to serve as source for disper- 
sal towards smaller islands [56,57]. 

The small population size of P. m. subniger as com- 
pared to that of P. m. manillae might be the result of the 
small size of the Polillo Islands combined with the loss of 
adequate habitat. This would explain the inferred recent 
bottleneck and the low values of genetic diversity of P. m. 
subniger (Table 5). Again, it should not be overlooked 
that our sample size for the subspecies is rather small. 

Aceros waldeni 

In the haplotype network (Figure 2C), all three haplo- 
types of the five individuals from Guimaras cluster 
among the Panay haplotypes. Unfortunately, we do not 
have microsatellite data for these museum specimens; 
however, mtDNA (Table 5) suggests a certain degree of 
genetic connectedness between these two islands. The 
water channel between Guimaras and Panay is nowadays 
about 2 km wide and obviously did not act as an insur- 
mountable barrier to gene flow. Between these two 
populations of A. waldeni, we did not find any major 
phylogeographic break. Nevertheless, each of the three 
haplotypes is unique. Thus, this island hosted independ- 
ent (yet closely related) mitochondrial variants that 
probably got lost when the population went extinct. 

However, despite having analyzed only a single mu- 
seum specimen from Negros (BMNH 1897.5.13.471), 
similar to what was observed in P. panini, this individual 
indicates a separation between the Negros and the 
Panay/Guimaras population (Figure 2C). The same pro- 
cesses described for P. panini above may thus also have 
had an impact on A. waldeni. The Negros population is 
believed to be extinct or nearly so, as no recent sightings 
have been reported [EC pers. obs.]. Thus, this differen- 
tiated mitochondrial lineage has already gone extinct or 
is at least threatened with extinction, if any remaining 
hornbills are still alive. 

Conclusions 

P. manillae is genetically more structured at the mtDNA 
level than P. panini. The latter is undergoing a very 



Sammler et al. BMC Evolutionary Biology 201 2, 1 2:203 
http://www.biomedcentral.com/1471-2148/12/203 



Page 13 of 14 



rapid continuing decline, associated with genetic deple- 
tion. Appreciable populations survived only on Panay 
and Negros. These two populations are genetically dif- 
ferentiated, as are those extinct or almost extinct from 
Guimaras and Masbate. This means, at least in the case 
of Guimaras, that part of the original genetic diversity of 
the species has already got lost. 

At the microsatellite level, we found evidences of a recent 
bottleneck for P. panini from Negros, but also for P. m. 
subniger (Polillo Islands). While the latter population might 
be small in size because the islands it colonizes are also 
small, Negros is a large island and its population shrank 
solely because of loss of adequate habitat. 

For A. leucocephalus, our data suggest that this species is 
(or was) genetically highly structured. A. waldeni has very 
rapidly decreased in number [12]. It survives with a sub- 
stantial population size only on Panay, i.e., in the forest of 
the Central Panay Mountain Range [23]. A. waldeni gener- 
ally exhibits lower genetic diversity values, both at microsa- 
tellites and mtDNA, than its less endangered sister species 
A. leucocephalus. Nevertheless, no sign for a recent bottle- 
neck was found for A. waldeni. For the extinct population 
on Guimaras we could show that an independent 
(yet closely related) lineage probably has already got lost. 
From what we found in P. panini,we are inclined to 
hypothesize that the extinct or almost extinct population of 
A. waldeni on Negros is most likely genetically differen- 
tiated from that still extant on Panay. This hypothesis is fur- 
ther supported by the mt haplotype of the single museum 
specimen from Negros. 

To sum up, the populations of the analyzed hornbill spe- 
cies are generally structured by the geological history of the 
area, although we found a few indications for across-island 
dispersal. By comparing patterns of genetic divergence and 
variability in both endangered Visayan hornbills species to 
that of their (less endangered) sister taxa, we revealed lower 
genetic diversity associated with their dramatic population 
decline and the extinction of genetically differentiated 
populations. Conservation efforts maintaining the surviving 
wild birds of virtually extinct populations are therefore par- 
ticularly important, as well as the preservation of their 
genetic potential in captivity. This applies especially to the 
potentially still existing populations of A. waldeni on 
Negros and P. panini on Masbate. 
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